Dynamic correlations in symmetric electron-electron and electron-hole bilayers 
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The ground-state behavior of the symmetric electron-electron and electron-hole bilayers is studied 
by including dynamic correlation effects within the quantum version of Singwi, Tosi, Land, and 
Sjolander (qSTLS) theory. The static pair-correlation functions, the local-field correction factors, 
and the ground-state energy are calculated over a wide range of carrier density and layer spacing. 
The possibility of a phase transition into a density-modulated ground state is also investigated. 
Results for both the electron-electron and electron-hole bilayers are compared with those of recent 
diffusion Monte Carlo (DMC) simulation studies. We find that the qSTLS results differ markedly 
from those of the conventional STLS approach and compare in the overall more favorably with the 
DMC predictions. An important result is that the qSTLS theory signals a phase transition from the 
liquid to the coupled Wigner crystal ground state, in both the electron-electron and electron-hole 
bilayers, below a critical density and in the close proximity of layers (d < ?'sfflo)j in qualitative 
agreement with the findings of the DMC simulations. 

PACS numbers: 71.10.-w, 73.21.-b,73.20.Qt 



I. INTRODUCTION 

In recent years, there has been considerable interest 
in the study of systems composed of two or more (equi- 
spaced) electron layers. The advances in the nanoscale 
semiconductor fabrication technology (such as the molec- 
ular beam epitaxy, the lithography techniques, etc.) have 
made available these electron systems in the coupled 
semiconductor quantum-well structures with a good con- 
trol on the electron number density and the interwell 
spacing. A variety of new interesting phenomena have 
been observed, due entirely to the presence of interlayer 
Coulomb interactions. The stability of new fractional 
quantum Hall statesEJ and tie discovery of an insulating 
Wigner crystal (WC) phasecl, in the bilayer electron sys- 
tem in the presence of a perpendicular magpetic field, 
are some prominent examples. TheorsticalD'El as well as 
recent diffusion Monte Carlo (DMC)|j'El studies have pre- 
dicted, even in the absence of magnetic field, the stabi- 
lization of the WC phase in the electron bilayer in the 
close proximity of layers d ^ Tga^ (see below for the defi- 
nition of relevant parameters) at sufficiently low electron 
density. More precisely, the critical density for the phase 
transition is predicted to shift towards the higher density 
side as compared to the corresponding value in the case 
of an isolated electron layer. In the following we shall be 
confining our discussion to the double layer system. 

Central in understanding the behavior of the layered 
electron systems are the intra- and interlayer many-body 
correlation effects. Since the early random-phase ap- 
proximation (RPA) study by Das Sarma and MadhukaiQ, 
there have been many investigations focusing on the role 
of many-body correlations. The theoretical techniques 
used have mostly relied on the extensions to two [dimen- 
sions (2D) of existiag theories. Zhang and Tzoaift, Neil- 
son and co-workersti, and Zheng and MacDonaldly have 



used the mean-field apprnximation of Singwi, Tosi, Land, 
and Sjolander (STLS)E2l to study the effect of correla- 
tions on the various ground-state properties of an elec- 
tron bilaye&i system. On the other hand, Kalman and 
co-workergl3 have incorporated correlations by satisfying 
the third-frequency sum rule of the density-density re- 
sponse function to study the collective modes. A general 
conclusion has emerged that the interlayer interactions 
add further to the importance of many-body correlations, 
which are already kaown to be very important in an iso- 
lated electron layeill3. Both in the STLS and Kalman's 
group approaches, correlations enter in the theory in the 
form of a static local-field correction (LFC). In STLS, 
the LFC's are obtained numerically in a self-consistent 
way, while in the latter approach their calculation rely 
on the knowledge of accurate pair-correlation functions. 
The extent of validity of the STLS or other theories in 
the present context can be tested by making a direct 
comparison with the accueate DMC simulation results of 
Rapisarda and SenatoreErQ. However, no such compari- 
son of the STLS results has been made so far. A compar- 
ative studji. is available only in case of an isolated elec- 
tron layerliJ, where it has been found that the STLS the- 
ory, although providing a significant improV|eHrent over 
the lower order random-phase and Hubbarcii3 approxi- 
mations, yet it fails to give an adequate description of 
correlations beyond > 3. In particular, it yields nega- 
tive value for the pair-correlation function (an unphysical 
result) at small separation for rg > 3. This failure of the 
STLS approach also appears in the bilayer problemEj. 
As usual, here Vg = 1 / {d^^^/rm) is the dimensionless den- 
sity parameter, with = h'^/{m*e'^) the effective Bohr 
atomic radius and n the in-layer areal number density. 
TO* is the effective (band) mass of electron. We recall 
that Ts also provides a rough estimate of the in-layer cou- 
pling, as ratio of the (independent particle) potential and 
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kinetic energies of the system. The failure of the STLS 
theory at higher Tg has been traced back to its basic as- 
sumption of treating as static (i.e., time- independent) the 
electronic exchange-correlation holell3. 

In the two-layer system, the interlayer interactions will 
result in an increase in Coulomb coupling as compared to 
that in an isolated layer at the same density. Therefore, 
we anticipate that the dynamics of correlations will be 
even more crucial in the two-layer system, with respect 
to an isolated layer. This forms part of our motivation for 
the present work. Here, we intend to examine the ground- 
state behavior of the symmetric electron-electron (e-e) 
bilayer by including dynamic electron correlation effects. 
To this end, we make use of the dynamic or quantum 
version of the STLS approach (qSTLS). The qSTLS the- 
ory embodies correlations beyond the conventional STLS 
approach and as an important improvement, its LFC is 
frequency-dependent. It should be pointed out here that 
the so-called qSTLS tlppry was originally developed by 
Hasegawa and ShimizutZl for the 3D electron system, and 
that its predictions of many-body properties, as adjudged 
by comparison with the MC results, are better irrespec- 
tive of the dimensionalityliJ and the carrier statisticsEJ. 
In view of the very recent DMC study of the symmetric 
electron-liole (e-h) bilayer by De Palo, Rapisarda, and 
Senatorell^, we also employ the qSTLS theory to exam- 
ine the ground-state properties of such system. For both 
the e-e and e-h bilayers we present results for the pair- 
correlation functions, the LFC factors, and the ground- 
state energy over a wide range of density and layer spac- 
ing. We also look for signs indicating an instability of 
the homogeneous liquid with respect to inhomogeneous 
phases of the charge-density-wave (CDW) and WC type. 
Finally, we compare our results with those of the DMC 
studies. 

The rest of the paper is organized as follows: In Sec. II, 
we present in brief the qSTLS formalism for the double- 
layer system. Results and discussion are presented in 
Sec. III. In Section IV, we conclude the paper with a 
brief summary. 



II. THEORETICAL FORMALISM 
A. Model 

We consider a double quantum-well structure with d 
as the center-to-center well separation. The carriers are 
electrons in one well and electrons or holes in the other, 
respectively for the e-e and e-h bilayer. The motion of 
carriers is free along the xy-plane and under the action 
of a double well potential profile in the z-direction. We 
assume that the wells are extremely narrow and the po- 
tential barriers along the z-axis are high enough so that 
the particles occupy only the lowest energy subband for 
the z motion and there is negligible overlap between the 
wave functions of particles in the two wells. The wells 
are assumed to be identical in each respect except for 



the charge of carriers in the e-h bilayer. Further, the 
bilayer system is assumed to be embedded in a uniform 
charge neutralizing background. On neglecting the effect 
of integrating over the finite extent of the particle wave 
function in the z-direction, the Coulomb interaction po- 
tential among the carriers is obtained as 



Vu'{q) = au'V{q)e"^^'~''^'', 



(1) 



with I = 1,2 the layer index and V{q) = 27re^/((7eo) the 
intralayer interaction potential. Above eo is the back- 
ground dielectric constant and = 1 and (— 
respectively, for the e-e and e-h bilayers. 

Apparently, the ground state of the above bilayer 
model will depend, apart from , on the interlayer spac- 
ing d. It turns out convenient to introduce an addi- 
tional coupling parameter as the ratio of the typical inter- 
layer and in-layer Coulomb energies, namely, 7 = rsOg/d. 
Thus, at T = 0, which is the case considered here, the 
bilayer model may be completely specified by and 7 
(or d). 



B. Density response function 

In the dielectric approach, the density-density linear 
response function x{q^^)^ which describes the response 
to an external potential V'^^*{q,Lo) that couples to the 
particle density, plays the role of a central quantity in de- 
termining the many-body properties of the system. For 
the bilayer, the linear response matrix ^/('(g, is for- 
mally defined by 



<5p,(9,^) = E^"'('?'^)^'''*('^''^)' (2) 
v=i 

where 5pi{q, uj) represents the induced particle density in 
the Zth layer and V{'''^'^{q,Lo) the potential that couples to 
the density in the layer /. For completeness, we give in the 
following a very brief account of the qSTLS formulation 
of the density response for the two-layer system. 

The equation of motion for the one-particle Wigner 
distribution function (WDF|)-|/;'^(r, p; t) (the superscript 
cr is the spin index) involvcsEil, the unknown two-particle 
WDF fii7 (r, p; r', p'; t). Progress can however be madfi 
resorting to the STLS approximate decoupling ansatzO, 
which in the two-layer system becomes 

/r' (r, p; r', p'; i) = /r(r, p; t)/,?' (r', p'; t)5;7' (|r - r'|), 

(3) 

with Qii? (|r — r'l) the equilibrium static pair-correlation 
function between carriers of spin a and a' in the layers 
I and V . Expressing the particle density in terms of the 
the one-particle WDFtZl, one readily obtains the induced 
density Spi{q,Lu) in the Ith layer as 

Spi{q,u;)=xnq,^) \vr'{q,u;) + Vr\q,u;)] , (4) 
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with the polarization potential, x'liQi'^) the density response 

function of nopiinteracting electrons in layer I (i.e., the 



Gu'{q,uj) 



1 /- rfq^ xl'(q,q';^)^r(gO 
I 



[5n'(|q-q'|)-&'], 



(6) 



the dynamic LFC factor that accounts for correlation ef- 
fects among carriers in the layers I and I'. In Eq. (|^), 
Su>{q) is a static structure factor and Xj'(q, q'; tj) the in- 
homogeneous Stern function given by 



x°(q, q';^) 



dk /0(k + q72)-/0(k-qV2) 



(2^)2 



w ■ 



?ik • q/m -|- ir; 



(7) 



where //"(k) is the usual non-interacting Fermi-Dirac dis- 
tribution function and ?7 is a positive infinitesimal. For 
q' = q, x?(qil';'^) reduces to Stern function 
Using Eqs. (p|), (^, and (||) the elements of the inverse 
of the linear response matrix are readily obtained as 



5w 



Vu'iq)[l-Gu,iq,uj)]. (8) 



ground-state energy theore: 
as 



Egs — Eq 



'2-'' to the two-layer system 



(11) 



where Eq = p'j^/{Aml) is the kinetic energy per particle 
of the non- interacting system (pp = hqp, is the Fermi 
momentum), A is the strength of Coulomb interaction, 
and i?*"*(A) is the interaction energy per particle given 
by 



1,1' 



dq 

(2^ 



XVu,iq)[Su>iq;X)-6iu]. (12) 



In the next section, we present results for the ground- 
state properties of the e-h and e-e bilayers. 



The fluctuation-dissipation theorem, which relates the 
static structure factors with the imaginary part of the 
linear response functions as 

Sii'{q) = / dujlmxu'{q,uj), (9) 

mr Jo 

closes the qSTLS set of equations for the density response 
matrix. Evidently, the response function calculation has 
to be carried out numerically in a self-consistent way. In 
view of the symmetry of the bilayer, we will have An — 
A22 and A12 = A21, where A refers to the general layer 
property. 



C. Pair-correlation function and ground-state 
energy 

The pair-correlation function gui (r) can be obtained 
directly from the inverse Fourier transform of the static 
structure factor as 

gii'ir) = 1 + ^ / - -^"'l- (10) 

The ground-state energy Eg^ (defined here as per parti- 
cle) is determined by a straightforward extension of the 



III. RESULTS AND DISCUSSION 
A. Pair-correlation functions 

Equations (|), (|), and (|) are solved numerically in 
a self-consistent way for Sw{q). The ^-integration in 
the computation of Sii'{q) (Eq. (^) is performed along 
the imaginary w-axis in order to avoid the problem of 
dealing with the plasmon pplf3s which anpear on the real 
w-axis (see de Freitas et alE3 and Refo). We accepted 
the solution when the convergence in Sii'{q) at each q 
in the grid of g-points was better than 0.001%. It is 
important to point out here that in all our calculations 
for the e-h bilayer we have taken m,* /m* = 1 ( m,* is the 
effective mass of hole) . 

Figures |^ ((a)-(g)) and || ((a)-(b)) show results for 
the intra- and interlayer pair-correlation functions, gii{r) 
and (712(7'), for the e-h and e-e bilayers, respectively, at 
Vg (< 10) and d values where the DMC results are avail- 
able for comparison. In order to have a close compar- 
ison between the qSTLS and STLS results, the STLS 
curves are also depicted in the same figures. We first 
discuss the results for the e-h bilayer: Looking at Fig. |l| 
((a)-(g)), we infer immediately that the pair-correlation 
functions in qSTLS are in overall better agreement with 
the DMC data than those in STLS. Among the notable 
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FIG. 1: (a)-(g) Pair-correlation functions g\\(r) (thin lines) 
and giiir) (thick lines) for the electron-hole bilayer at differ- 
ent values of Ts and I/7 = d/rsoZ- Solid and dashed lines are, 
respectively, the oS'HjS and STLS results; the prediction of 
DMC simulationsLJO for 311 (r) (o) and (712 (f") (+) are also 
shown, (h) i7ii(r-) (thin lines) and gi2{r) (thick lines) at in- 
dicated rs and I/7 values, according to qSTLS; STLS results 
are shown by dash-dot lines. 



features, the qSTLS theory, as a marked improvement 
over STLS, accounts fairly well for the oscillatory struc- 
ture that develops in the DMC gii{r) and gi2{r) with 
increasing rs, both in terms of amplitude and period. 
The qSTLS 511 (r) satisfies the positive definiteness cri- 
teria of probability for rs up to 10, whereas the STLS 
9ii{f) becomes slightly negative at small r for rs > 3. 
We also notice that, common at each rs, the quality of 
agreement between theory and DMC data somewhat di- 
minishes for increasing values of 7. Moreover, quite sim- 
ilar is the trend of agreement with increasing rs at a 
given value of 7; for instance, compare 511 (r) and .912 (f") 



FIG. 2: (a)-(b) Pair-correlation functions (?ii(r) and (712 (r) 
for the electron-electron bilayer; Curves are labeled as in Fig. 
|l} The results of DMC simulations are from Refs. 



at 1/7 — 0.5 and Ts = 2,5,10. This shortcoming of 
the qSTLS theory seems to stem from the fact that the 
frequency-dependence in its LFC's represents a quantum- 
mechanical correction to the conventional STLS theory, 
while the dynamics of spatial correlations among carri- 
ers, which is expected to become vital at higher values 
of Coulomb coupling, is still missing. It can in fact be 
showntJ that the qSTLS LFC's reduce formally to the 
frequency-independent STLS LFC's in the limit fi ^ 0. 

Since decreasing d at a given Ts or increasing Ts at a 
given d result in an increase in the Coulomb coupling 
among carriers in the bilayer, the qSTLS assumption of 
static spatial correlations (i.e., the assumption of using 
static pair-correlation function in Eq. (H)) is expected 
to become relatively less reliable at larger values of 7 or 
Ts- Further, we find that it becomes almost impossible 
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TABLE I: The ground-state energy Egs per particle (in units 
of effective Rydberg) at different and I/7 values (7 — 
VsaQ/d) for the e-h bilayer, according to qSTLS and STLS. 



rs 1/7 


qSTLS 


DMC 


STLS 


2 0.1 


-0.6519 


-0.6947 




2 0.2 


-0.5831 


-0.6116 


-0.5757 


2 0.5 


-0.5412 


-0.5405 


-0.5307 


5 0.5 


-0.3057 


-0.3125 


-0.3010 


K in 
l.U 


-0.3015 


-0.3009 


-0.2970 


10 0.5 


-0.1724 


-0.1801 


-0.1706 


10 1.0 


-0.1697 


-0.1715 


-0.1682 


10 1.5 


-0.1695 


-0.1703 


-0.1680 


TABLE II: Theg 


;round-state energy 


Egs per particle (in units 


of effective Rydb 


erg) at — 10 for different I/7 values (7 — 


rsag/d) for the c 


i-e bilayer, according to qSTLS and STLS. 


DMC results are 


from Ref. § 






rs 1/7 


qSTLS 


DMC 


STLS 


10 0.5 


-0.1707 


-0.1781 


-0.1699 


10 1.0 


-0.1685 


-0.1713 


-0.1677 


10 1.5 


-0.1683 


-0.1705 


-0.1675 



to obtain the self-consistent solution in both the qSTLS 
and STLS above a critical value of 7 (r^) at a given 
(7). The critical parameters are of course different in 
the two approaches. This is the reason that the STLS 
curves are absent at rg = 2 for I/7 = 0.1 in Fig. |^ (c). 
The difficulty in obtaining the self-consistent solution, as 
we will see in detail in a subsequent section, is related 
to the instability of the system against transition to a 
density-modulated ground state. 

For the e-e bilayer, the DMC pair-correlation functions 
are at present available only in the strong coupling region 
fs > 10. Figure ^ ((a)-(b)) present a comparison of our 
results at = 10 for I/7 = 1 and 0.5. We notice that 



the qSTLS provides a reasonable estimate for 511 (r), but 
it fails ( together with STLS) to give a satisfactory de- 
scription of gi2{r) at 1/7 = 0.5. This points once again to 
the importance of dynamics of spatial correlations among 
electrons in the e-e bilayer. The improvement over the 
STLS predictions is again quite noticeable. 



B. Ground-state energy 

The self-consistently obtained static structure factors 
5ii(g; A) and S'i2(g; A) are used in Eq. ( p^ ) to calculate 
E™*-{\) as a function of A. The ground-state energy is 
then determined by performing the coupling-constant (A) 
integration in Eq. (pT|). Results for the ground-state 
energy per particle for the e-h and e-e bilayers are given, 
respectively, in Tables | and || at (< 10) and 7 values 
where DMC data are available for comparison. The STLS 
results are also reported. Apparently, the qSTLS results 
compare more favorably with the DMC data. There is 
an increase in error with respect to the DMC data with 
increasing Vg at a given 7 and with increasing 7 at a given 
Ts , which obviously is the reflection of the behavior of the 
pair-correlation functions under these conditions. 



C. Density-modulated ground states 

The DMC studies have predicted that both the e-e 
and e-h bilayers will favor energetically the WC ground 
state above a critical value of the in-layer coupling, 
at given d. We calculate here the static {uj — 0) gen- 
eralized susceptibility (i.e., density-density response) in 
the liquid phase to find out any evidence that the qSTLS 
theory might provide for the transition to a WC ground 
state. If such a transition does occur, it may appear in 
the static susceptibility as a divergence at the reciprocal 
lattice vector (RLV) of the WC lattice. Diagonalizing 
the density response matrix (H), the static susceptibility 
is obtained as 



X±(9,0) = 



X?(9,0) 



1-X?(g,0) [Fii(g)(l 



Gii(g,0))±Fi2((z)(l 
I 



Gi2(g,0))]- 



(13) 



The -|- and — signs correspond, respectively, to the in- 
phase and out-of-phase (tt) modes of density modula- 
tions 5p{q, 0) in the two layers. An inspection of Eq. 
( |l3| ) makes it clear that x±(<7, 0) can exhibit divergence 
at some q value only if the quantity within the square 
brackets in its denominator becomes sufficiently nega- 
tive (on account of the negative sign of Xi(g,0)) . In 
an isolated layer (i.e., when Vi2{q) — 0) this can happen 



only if GaUci, 0) has values exceeding unity. However, in 
a bilayeiuij, the interlayer interaction term can cause a 
divergence even if Gii((;, 0) has values below unity. Ap- 
parently, it is the in-phase component of susceptibility 
that can have divergence in the e-h bilayer, while it is 
the out-of-phase component in the e-e bilayer. Further, 
as the e-h and e-e correlations are of opposite nature (i.e., 
Gi2{q, 0) is negative in the e-h bilayer, while it is positive 



7 



in the e-e bilayer), they will act, respectively, to support 
and oppose the formation of density-modulated phase (if 
any) in the e-h and e-e bilayers. 

The static LFC factors, Gii(g,0) and Gi2(<z, 0), re- 
quired in the calculation of the susceptibility are de- 
termined by using the self-consistently obtained static 
structure factors in Eq. (||). We find quite generally 
that x+{Qt^) (X-(9:0)) exhibits for the c-h bilayer (e- 
e bilayer) a strong peak-structure at a finite wavevector 
value in the close proximity of two layers {d ^ r^aQ). A 
critical layer separation dc (critical r^) at a given (d) 
is encountered below (above) which it becomes almost 
impossible to obtain the self-consistent solution. Trac- 
ing carefully the different steps involved in the solution 
of the qSTLS equations, the difficulty in obtaining the 
convergent solution is found to be related directly to the 
emergence of the strong peak structure in x±((Z, 0). For 
d < dc Sit a, given or for > at given d, a numeri- 
cal instability (singularity) appears in Sw (q) during the 
iterative calculation, at a g coinciding exactly with the 
peak position in Y±(g,0), while calculating Siii{q) from 
Xw {q, w) in Eq. (g). We find it extremely difficult to han- 
dle the instability inJjhe self-consistent calculation and 
therefore, are unableEj to find the convergent Su> (q) be- 
low (above) a critical value of d (r^). This forbids us, in 
turn, to calculate x±iQT^) below (above) a critical value 
of d {ts) at a given Ts (d). On the other hand, the emer- 
gence of a strong peak in x±{Qj 0) a-t finite wavevector 
value followed by the numerical instability at the same 
wavevector during the self-consistent calculation of the 
density response functions could be interpreted in our 
theory as an indication for the onset of a phase tran- 
sition from the liquid to the density-modulated ground 
state. 
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1. Electron-hole bilayer 

Results for x+(<Z, 0) at some selected values of and 
d are shown in Fig. || ((a)-(c)). For < 5, x+{QiO) 
has a single strong peak at small q; for instance, at 
rg = 2 (Fig. 11(a)) the peak is positioned at q/qp ~ 0.2 
and dc/irsCLQ) ~ 0.095. For > 5, however, a second 
peak starts developing in x+{l-i^) a-t q/qr ~ 2.5, with 
its strength relative to that of the peak at small q/qp 
growing continuously as a function of increasing r^. At 
Ts = 10 (Fig. ||(b)), the peak at q/qp ~ 2.5 eventu- 
ally dominates the small-g peak and dcKvsa'^) ~ 0.36. 
With further increase of r^, the small-g peak disap- 
pears completely as shown in Fig. |^(c) at rg = 15; 
dc/irsO-o) ~ 1.10. The position of the peak located at 
q/qp Ri 2.5 matches quite closely with the RLV for a 
triangular Wigner lattice (« 2.7) and therefore, we spec- 
ulate that this peak signals a transition to the coupled 
in-phase WC ground state in the e-h bilayer. Our spec- 
ulation draws further support from the fact that near 
the transition point gii{r) and 5i2('') exhibit strong in- 
phase oscillations typical of an ordered phase, and this 
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FIG. 3: (a)-(c) In-phase component of the static density sus- 
ceptibility x+iljO) for the electron-hole bilayer at different 
values of ra and I/7 — d/raag, according to qSTLS. Legends 
indicate the values of I/7; STLS results are shown for com- 
parison by dash-dot lines in (c) at rs = 15 and indicated I/7 
values. 



feature of g{r) is illustrated in Fig. |l|(h) at = 15. 
On the other hand, the small-g peak, whose position 
varies with r^, indicates the instability of the e-h liq- 
uid against a CDW ground state. Thus, there seems to 
occur a crossover from a CDW to a WC ground state at a 
critical density « 10. The qSTLS indication of a tran- 
sition to a WC state is in agreement with the findings of 
the DMC studyO. However, the critical Vs for the onset 
of the phase transition is underestimated by a factor of 
about 2. The qSTLS theory signals a Wigner crystal- 
lization also for an isolated electron layer, at « 17, 
again with an overestimate of the critical coupling by a 
factor of abjMit|22i as compared with the available DMC 
predictionsElBE^O. Nevertheless, an important result to 
note is that the electron-hole correlations in the e-h bi- 
layer act to lower the critical value, as compared to 
that in an isolated layer, by a factor of about 1.7, which 
matches closely with the DMC prediction. 

Furthermore, the DMC studyllj has found that the WC 
ground state (also, the liquid state when < r^) be- 
comes energetically unstable against transition to an ex- 
citonic ground state as the layers of electrons and holes 
are brought close to each other (d < TsCIq)- There is some 
indication in our results, in terms of the steady buildup 
of .912 — 0) in tjae close vicinity of layers, for the for- 
mation of excitonaHj, but there is no apparent way in our 
theory to directly detect the transition to the excitonic 
ground state. 

It is appropriate at this point to draw a comparison of 
our results with the previous werk on the e-h bilayer by 
Liu et aB and Szymanski et al3. Liu et al treated the 
correlations within the STLS approach and found at all 
Tg, in the close vicinity of layers, always an instability 
towards the CDW ground state. The STLS x+((?,0) is 
plotted for comparison at = 15 in Fig. |3|(c); the CDW 
wavevector q/qp is « 2 and dc/ (j'soj) ~ 0.48. Also shown 
for comparison in Fig. ^h) are the STLS gii{r) and 
5i2(r) near the CDW instability (at d/(rsal) = 0.49) at 
rs = 15. On the other hand, Szymanski et al employed 
the STLS local-fields to include the correlation effects, 
but instead of carrying out the fully self-consistent STLS 
calculation, they fixed the intralayer LFC GiJ,n) through 
the MC structure factor for an isolated layered, while the 
interlayer LFC G'i2('z) was determined self-consistently 
by keeping Gii(q) as a fixed input. Following this mixed 
STLS procedure, they predicted the presence of both the 
CDW and WC instabilities, and also a crossover from the 
CDW state to the WC state at rs « 15. Thus, there is a 
qualitative similarity between the qSTLS results and the 
findings of Szymanski et al. 

It is gratifying to note that the qSTLS theory seems 
to capture at least qualitatively the WC transition in 
the e-h bilayer. The breakdown of the STLS approach 
even at the qualitative level seems to have roots in its as- 
sumption of treating correlations through the frequency- 
independent (i.e., static) LFC's. To elucidate this view- 
point, the static (w = 0) qSTLS LFC's and the STLS 
LFC's, which in fact determine completely the behav- 
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FIG. 4; Intralayer (in panel (a)) and interlayer (in panel (b)) 
static local-field corrections G'ii(g, 0) and Gi2(g, 0) for the 
electron-hole bilayer at rs = 15 and different layer spacings, 
according to qSTLS. Legends indicate the values of I/7 = 
d/rsOo; STLS local-fields are also shown. 



ior of x+('7, 0), are compared in Fig. |j at rs = 15 for 
different d's. The LFC's differ markedly in the two ap- 
proaches. In particular, we note that Gii{q,0), in con- 
trast with the STLS Gii(g), exhibits an oscillatory be- 
havior, with pronounced maximum at q/qp ~ 3 and 
with its values lying well above unity in the relevant 
wavevector region of 2 < q/qp < 3, before saturating 
to its limiting value of (1 — 511(0)). Also, Gi2{q,0) is 
rapidly varying in the same g— region. Consequently, 
the qSTLS effective static intralayer interaction, namely 
[Vii(g)(l — Gii{q,Oj)], becomes attractive in this q- 
region, and this effect of Gii(g, 0) in combination with 
the attractive interlayer e-h correlations gives rise to a 
strong peak in x-i-((7, 0) at q/qp ~ 2.5. The static LFC 
factor exhibits a similar oscillatory behavior in an iso- 
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FIG. 5: (a)-(b) Comparison of the ground-state energy per 
particle Egs between the unpolarized and polarized phases of 
the electron-hole bilayer at d/aS = 3.6 (in panel (a)) and 20 
(in panel (b)), according to qSTLS. 



lated electron layer. Though the peaked-structure in the 
qSTLS static LFC factor seems to enable the qSTLS the- 
ory to capture quahtatively the WC instability, we have 
to mention that at least for an isolated layer such an 
oscillatory structure is exaggerated as compared to that 
found in QMC simulationg23. 

We also examine the stability of the ground state of 
the e-h bilayer as a function of spin-polarization. For 
simplicity reasons, only two states of spin-polarization, 
namely the fully polarized and unpolarized (which oth- 
erwise is the case throughout the paper), are considered. 
A comparison between the ground-state energies of the 
e-h bilayer in the two states of spin-polarization (Fig. ||) 
reveals that a polarization transition occurs at a critical 
density in the liquid phase from the unpolarized to the 



FIG. 6: Critical layer spacing dc/ (tsOq) (at a given Va) for the 
transition from the liquid to a density-modulated phase as a 
function of in the unpolarized (triangles A) and polarized 
(circles o) phases of the electron-hole bilayer, according to 
qSTLS. For each case the arrows show the critical density 
where crossover from the charge- density- wave instability to 
the Wigner crystal instability occurs; The lines are just a 
guide for the eye. 



polarized state before the unpolarized liquid could actu- 
ally make transition to the seemingly WC phase. The 
critical density for polarization transition depends only 
weakly on d; for instance the critical Ts decreases from 
6 to 5.5 as d/a^ increases from 3.6 to 20. We empha- 
size here that d/ap = 3.6 is the critical value of layer 
spacing at = 10 for the WC instability. In this per- 
spective, it becomes interesting to investigate the ground 
state of the polarized e-h bilayer. We find after analyz- 
ing the static susceptibility results that the qualitative 
behavior of the e-h ground state does not depend upon 
the spin-polarization. However, the crossover from the 
CDW phase to the WC phase now occurs at « 27. 
The dependence of the point of instability on the spin- 
polarization is depicted in Fig. |^. Evidently, the critical 
spacing in the polarized phase lies always below to 
that in the unpolarized phase. 



2. Electron- electron bilayer 

In contrast with the e-h bilayer, the interlayer correla- 
tions in the e-e bilayer tend to oppose the transition to 
the density-modulated phase. But, we find that this ten- 
dency of the interlayer correlations depends crucially on 
the rate of their growth with decreasing d, and that this 
rate is not strong enough to preclude transition to the 
density-modulated phase. We find indication for both 
the CDW and WC instabilities. 

The X- ily 0) results are shown in Fig. |^ at = 10 
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FIG. 7: (a)-(b) Out-of-phase component of the static density 
susceptibility X-[q, 0) for the electron-electron bilayer at = 
10 and 16 and different values of I/7 = d/raa*), according to 
qSTLS; legends give the I/7 values. 



and 16 for some selected d's. The CDW instability now 
completely dominates for up to 15. But, we notice an 
interesting behavior of X-(9i0) at rg = 16. X-('J'iO) ex- 
hibits a single strong peak at q/qp ~ 2.5 when the layers 
are widely separated. The peak-height initially increases 
with decrease in d, nearly diverges at d/{rsa'^) ~ 1.56, 
and then decreases monotonically with further decrease 
in d. There starts developing, however, a second peak 
in x-(9,0) at q/qp ~ 0.5 for d/{r sUq) < 0.5, and this 
peak eventually dominates at d/lrsa^) « 0.45 and then 
appears to diverge for d/lrga^) < 0.42. This might pos- 
sibly imply that there is a crossover to the WC ground 
state at « 16, with the WC state however remaining 
stable only over a certain range of layer spacings, with a 
transition back to the CDW-like phase when the distance 



FIG. 8: Pair-correlation functions gii{r) (thin lines) and 
gi2{r) (thick lines) in the electron-electron bilayer at = 16 
and different values of I/7 = d/rsa^, according to qSTLS. 
Legends give the I/7 values. Note that the gii{r) at 1/7 = 
1.88 and at I/7 = 1.56 are practically indistinguishable. 



is further diminished. The pair-correlation functions in 
the relevant range of d values are shown in Fig. ^. Evi- 
dently, there is a continuous decline in the amplitude of 
oscillation in 511 (r) with decreasing d/(rsa5) below 1.56 
pointing to the obvious screening of in-layer correlations 
by the interlayer interactions. The present suggestion 
of the stability of the WC ground state in the e-e bi- 
layer only at intermediate distawjes is compatible with 
the findings of DMC calculationsBEl, which for Ts not too 
large predict the WC phase to remain stable only at in- 
termediate layer spacing, with the liquid phase becoming 
stable again at smaller spacing. 

The STLS approach again does not give any indica- 
tion of the WC phase transition. As in the case of e-h 
bilayer, this inability of STLS is a manifestation of its 
assumption of frequency-independent LFC's. This point 
is illustrated in Fig. ^ by drawing a comparison between 
the static qSTLS LFC's and the STLS LFC's. We notice 
that Gii{q, 0) and Gi2{q, 0) both exhibit, in contrast with 
their STLS counterparts, a pronounced peaked-structure 
in the intermediate wavevector region before converging 
to their respective limiting values, which are close to the 
corresponding STLS limiting values. 



D. Dynamic local-fields 

We have seen above that it is the frequency-dependence 
of LFC's which brings in a marked difference in the 
qSTLS description of correlations as compared to the 
STLS one. Therefore, we examine the LFC's for their 
dependence on frequency. The results for the e-h and e-e 
bilayers are reported, respectively, in Figs. |l^ and |ll| at 
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FIG. 9; Intralayer (in panel (a)) and interlayer (in panel 
(b)) static local-field corrections Gii{q,0) and Gi2{q,0) for 
the electron-electron bilayer at = 16 and different layer 
spacings, according to qSTLS. Legends indicate the values of 
1/7 = d/rsOg; STLS local-fields are also shown. 



FIG. 10: Frequency-dependence of intralayer (in panel 
(a)) and interlayer (in panel (b)) local-fields Gii{q,uj) and 
Gi2{q,iij) for the electron-hole bilayer at q/qr = 1-5, Vs = 5, 
and different values of I/7 = d/rgaQ. Thin and thick lines rep- 
resent, respectively, the respective real and imaginary parts; 
legends indicate the I/7 values and S^f is the Fermi frequency. 



<lllF = 1-5, Ts = 5, and different d! s. The real and imag- 
inary parts of both the intra- and interlayer LFC's are os- 
cillatory functions of frequency - a_£eatuEe which is_analo- 
gous to that has been found in 3DEI, 2jM, and IDEI elec- 
tron systems. In the large frequency limit, both Gn (q, lo) 
and Gi2(g,w) approach formally the STLS LFC's. The 
LFC's in the e-h and e-e bilayers have qualitatively a sim- 
ilar dependence on frequency. However, we notice that, 
apart from the obvious difference of sign between the in- 
terlayer LFC's of the e-h and e-e layers, the attractive in- 
terlayer correlations in the former are relatively stronger 
in the close proximity of layers. Similar is the behavior 
of LFC's at other q and values. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have presented a study of the ground- 
state behavior of the symmetric e-e and e-h bilayers by 
including the effect of dynamic correlations within the 
qSTLS theory. We have found that the inclusion of the 
dynamical nature of correlations introduces quantitative 
as well as qualitative differences in the description of 
many-body properties as compared to static mean-field 
theories of the STLS type. The qSTLS predictions for 
the intra- and interlayer pair-correlation functions and 
the ground-state energy are found to be in overall bet- 
ter agreement with the DMC results. The growing os- 
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FIG. 11: (a)-(b) Dynamic local-fields for the electron-electron 
bilayer; Description of curves is exactly the same as in Fig. 
0. 



dilatory trends in the DMC intralayer correlation func- 
tion with increasing Ts are accurately reproduced both in 
terms of amplitude and period for Vg < 10 and I/7 > 0.5 
- a feature that is missing in the STLS results. How- 
ever, the degree of agreement with the DMC data, specif- 
ically at small interparticle separation, becomes some- 
what worse with increasing 7 = rsap/d, at a given or 
d, i.e., in the strong coupling regime. Another unique and 
important feature of the qSTLS theory is that, in both 
the e-h and e-e bilayers, it exhibits an instability towards 
a coupled WC ground state, below a critical density (i.e., 
for Vs > r^) and in the close proximity of the layers. 
Moreover, at high density (for rg < r^) it indicates tran- 
sition to a CDW ground state. Thus, a crossover from 
the CDW instability to the WC instability takes place at 
= r^. Our prediction of Wigner crystallization agrees 
qualitatively with the findings of the DMC calculations. 
However, the critical value of (i.e., r^) is underesti- 
mated. This discrepancy in the estimate of seems to 
reflect the misrepresentation of short-range correlations 
at large 7 or in the qSTLS theory. Since the Coulomb 
coupling among carriers grows with increasing 7 or and 
the exchange correlations are relatively less important in 
the strong coupling regime, we believe that the failure 
of the qSTLS theory in such a regime arises from its ne- 
glect of the dynamics of the Coulomb correlations, i.e., 
of the Coulomb correlation hole. Another contribution 
to the failure is possibly due to the neglect of interaction 
effects on the momentum distribution function. These 
issues deserve further investigation. 

Finally, we have also employed the qSTLS theory to 
examine the spin-polarization effects in the e-h bilayer. 
Interestingly enough, a polarization transition is found to 
take place from the unpolarized to the polarized liquid (at 

« 6) well before the unpolarized liquid could actually 
make transition to the WC ground state. The polarized 
e-h bilayer too supports the CDW and WC instabilities, 
but the crossover density is now lowered to rg ~ 27. 
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